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Backbone-induced semiconducting behavior in short DNA wires 
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We propose a model hamiltonian for describing charge transport through short homogeneous 
double stranded DNA molecules. We show that the hybridization of the overlapping tt orbitals in 
the base-pair stack coupled to the backbone is sufficient to predict the existence of a gap in the 
nonequilibrium current-voltage characteristics with a minimal number of parameters. Our results 
are in a good agreement with the recent finding of semiconducting behavior in short poly(G)-poly(C) 
DNA oligomers. In particular, our model provides a correct description of the molecular resonances 
which determine the quasi-linear part of the current out of the gap region. 
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The attempt to understand the mechanism of electron 
motion along DNA is the source of an intense debate in 
the biochemical and chemical physics communities [Q. 
Solving this problem is an essential step for the devel- 
opment of DNA-based molecular electronics. New in- 
sights to this issue are brought by recent breakthroughs 
in direct measurements through DNA molecules (|, [|, 
[|, ||, ^, 0, Transport measurements through nanos- 
tructured systems are potentially capable of addressing 
the basic issues of the conduction properties of molecular 
and supramolecular aggregates. The aftermath for the 
realization of molecular electronics devices is straightfor- 
ward H . It is thus not surprising that DNA molecules 
became the subject of an intense study concerning their 
potency to carry an electric current [g, ||, Q, ||, ^ [?], 
and to provide a scaffold for the metal assembling of 
highly conductive nanowires ^| . From a nanoelectron- 
ics perspective, the DNA possesses ideal structural and 
molecular-recognition properties, and the understanding 
of the charge transport through DNA may result in the 
ambitious goal of self assembling nanodevices with a def- 
inite molecular architecture ||Io|| . 

The hypothesis that double stranded DNA supports 
charge transport as a linear chain of overlapping tt or- 
bitals located on the stacked base-pairs, already ad- 
vanced in the early sixties Jll[ , received first experimental 
boosts only recently via long range electron transfer mea- 
surements . As far as transport through DNA is con- 
cerned, the available experiments arc still controversial 
mainly due to the complexity of the environment and the 
molecule itself (sequence variability thermal vibra- 
tions...). Concerning theory, the most reliable procedure 
to tackle these systems would be the ab initio quantum 
chemistry approach. However, massive numerical costs 
complicate its use for realistic biological systems Q . To 
our knowledge, at the present time, only few density- 
functional-theory (DFT) calculations for DNA molecules 
are available ^ [lf|. In a parallel development partic- 
ular aspects of the DNA transport phenomenology have 



been explained as mediated by polarons , solitons jl7) , 
electrons or holes Jl], . Such lack of a unifying theo- 
retical scheme calls for reproducible and unambiguous 
experimental results which are still a great technological 
challenge. 

Recently, Porath et al. [Q have reported nonlinear 
transport measurements on 10.4 nm long polyguaninc- 
polycytosine DNA, corresponding to 30 consecutive GC 
base-pairs, attached to platinum leads (GC-device). The 
measured room temperature current-voltage (I-V) char- 
acteristics show typical semiconducting features with a 
gap of the order of 1 V. Furthermore, the poly(G)- 
poly(C) DNA molecule has typical electronic features 
of a periodic chain, as the first DFT calculations have 
indicated Qj. This may support the idea that, differ- 
ently from natural A-DNA (complex sequence), where 
the sequence variability or the attachment to the surface 
could lead to electron localization over very few base- 
pairs ]l9| ], in short suspended GC-devices band-like con- 
duction might be the relevant transport mechanism. 

Motivated by such considerations, in this Letter, we 
introduce a minimal model for charge transport through 
GC-devices and show that the semiconducting behavior 
of the observed low temperature I-V curves can be ex- 
plained by the hybridization of the G-G tt stack with 
the transversal backbone reservoirs. The HOMO-LUMO 
structure, as estimated by Refs. |?], [2l], suggests that 
hole injection into the GC-devices might be fairly de- 
scribed using a tight-binding model by a hamiltonian 
comprising three terms H = H mo \ + H lcs ,ds + H coup \. 
We describe here a short poly(G)-poly(C) DNA molecule 
(N = 30 base-pair long) [|| as three-band model given by 

H mol = e^blb^-tu h l h 3° 

i,a,a— ± z,cr,a— ± 

where b\ a is the creation operators for charges with spin 
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Here A„ = — ImE;, is the spectral density of the metal 
molecule coupling. For notational convenience, we write 
down the relations for identical metal pads, El = Er = 
E. Gin is the molecular Green function between the two 
contact sites dressed by the lead self-energy. The cal- 
culation of Gin can be pursued analytically Q leading 
to 
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FIG. 1: Schematic view (left) of a fragment of poly(G)- 
poly(C) DNA molecule; each GC base-pair is attached to 
sugar and phosphate groups forming the molecule backbone. 
On the right side, the diagram of the lattice adopted in build- 
ing our model, with the n stack connected to the isolated 
states denoted as ±-edges. 



a in the G base site i = 1, . . . , TV), and c\ aa the one 
in the a-edge. The latter accounts for the upper sugar 
group sites and, possibly, for the C bases with the relative 
lower strand sites, as schematized in Fig. |l|. 

Charges can propagate along the n orbital stack via 
the nearest neighbour hopping probability in, or be hy- 
bridized to a-edges by £j_ Q ; Sb {s a ) is the energy level 
of the localized b (c a ) charge. We have fixed the num- 
ber of transversal paths per central site to two (a = ±). 
Models with more transversal hoppings are all equivalent 
(via a canonical transformation) . The numerical value of 
t± is then consequently renormalized by the number of 
such transversal paths. The leads are described by the 
hamiltonian 



,N+1 



N+l 



leads 



= E 



(1) 



k.o-,i/=L,R 



where k denotes the wave vector, and describes the 
single-electron dispersion relation of the i^-lead, measured 
with respect to the Fermi energy [i v (^l,r = ±eV/2). 
The coupling between the leads and the molecule ends 
can be described by a tunneling amplitude C/k of the elec- 
tron in the state (k, a, v) of the ^-lead to the molecule-end 
sites: 



H, 



coupl 



a Lh b la 



a kcrR^Wcr 



h.C. 



(2) 



The effect of the metal pad on the molecule is given by 
the self-energy 

E„ = 2^|C/ k | 2 G„ (k,E), 

k 

where G„ (k, E) is the retarded Green function for the 
isolated i^-lead, and the factor two accounts for the spin 
degeneracy. 

The transmission function, T = 4AlAr |<7iat| 2 , is ob- 
tained by making use of the Fisher-Lee relation p2~l. 



with £ 7V ($) = ($ + v 7 * 2 - iy~ ri - ($ - V$ 2 - 1) 

The relevant one-particle Green function at every site 
in the 7r stack has been renormalized according to the 
hybridization with the backbone states so that 



2*n 




(3) 



where Q v = (E + iQ + — e n ) are the bare (isolated) 
Green function for the three site classes (ry = b, ±). 

The backbone coupling, Eq. (||), controls the opening 
of a gap in the transmission. This can be intuitively un- 
derstood within the standard treatment for the leads. In 
fact, for bulky electrodes, in the wide band limit, ex- 
pressing the spectral density in units of tu, E = —i6t\\, 
the transmission can be written as 



T = 



4S 2 sin 2 (i9) 



(sin(/V + l)i9 - S 2 sin(JV - l)#y + 45 2 sin 2 iVt? 

T is an even function of 4> =: costf with N resonances 
at values of the molecular orbitals broadened by the di- 
mensionless lead spectral density 6. For simplicity, we 
consider here equal strengths in the backbone coupling, 
t± := tx+ = ^_L-; which does not imply any physical 
assumption if e v = for ?y = b (preserving the charge 
neutrality) and for r] = ± (no gating). In the absence 
of such a backbone coupling, $ reduces to the energy 
of the incoming charge relative to the G base on-site 
energy and in units of the band width of the n band, 
<5(i_i_ = 0) = (E — £;,)/2£||; the transmission of a iV-atom 
molecular wire is recovered ^3|. When a finite back- 
bone coupling is considered, the energy is renormalized 
through Eq. (|), thus a gap, At, is opened in the trans- 
mission following the lines of Fig. ||. In the absence of 
gating, the transmission gap reads 



A T = 2ji 



2t 2 



It 



(4) 



while the width of each of the two side bands is 2t 



and thus independent of tj_. This behavior can be 
also understood by referring to the dispersion relation, 
$(£■) = — cosg, in the limit of an infinite wire (N 3> 1), 
q being the longitudinal momentum in units of lattice 
spacing. The absence of electronic states between the 
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FIG. 2: Energy renormalization due to the backbone cou- 
pling. Different lines correspond to different values of t±/2t«, 
here e+ — e_ = et- The inset shows the case t± = 2t\\ 
where the two backbone are further gated with opposite sign 
e± = Eb ± tii. All curves are antisymmetric in <fr(t± = 0) = 
(E — eb)/2t\\. 

two emerging branches, the HOMO-LUMO gap of our 
model hamiltonian, determines the gap in the transmis- 
sion probability. As to strengthen the intuition of this 
high reflectivity near zero energies, one can regard our 
model as the extreme discretization of a phase coherent 
quantum waveguide with a fish bone shape (right hand 
side of Fig. @), where a low energy incoming particle has 
a high probability to be localized on the side transepts. 

The calculation of the current can be pursued within 
the scattering formalism j24[ 

I=^J dET(E) (f h (E) - ME)) , (5) 

where the Fermi functions f v (E) = l/(e! 3 ( E ~ ,i,u ' + 1) are 
controlled by the lead electrochemical potentials /iL,R- 
Eq. (||) is a reasonable estimation for the truly nonlin- 
ear current when the bridge system is a finite molecular 
chain. This has been recently shown by noncquilibrium 
Green function calculations (Keldysh formalism) through 
oncdimensional dot arrays [E5|. Moreover, since we con- 
sider here the low temperature current-voltage measure- 
ments from the experiment in Rcf. || (two representative 
examples are plotted in Fig. ||), the fit to the current can 
be tuned by comparing the experimental differential con- 
ductance with the transmission function (the upper inset 
of Fig. |^ shows the calculated transmission for the data 
displayed in blue). 

Theoretical results, plotted as solid lines in Fig. ||[ show 
a good overall agreement with the description of the gap 
and the molecular energy levels along the almost linear 
part of the experimental I-V characteristics. Theoreti- 
cal parameters have been obtained by a \ 2 minimization 
over i|| , constraining t± to give the observed experimental 
current gap. The latter is principally induced by the cou- 
pling t± to the side sites (see Eq. (0)) and reproduced into 
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FIG. 3: Low temperature I-V characteristics of two typical 
measurements at 18 K (blue circles) and at 3.6 K (red circles). 
Solid lines show the theory curves following the experimental 
data. The insets show the transmission calculated after the 
blue data (upper) and the normalized differential conductance 
(lower). The parameters used are ty = 0.37 eV and t±± = 
0.74 eV for the blue measurement, and t\\ = 0.15 eV and 
t±± = 0.24 eV for the red one §§. 

a gap in the current-voltage curves as a result of the inte- 
gration in Eq. (^J). For a fixed value of t\\ the current gap 
is an increasing function of t±. That is why the smaller 
hopping parameters used for the fit to the red data in 
Fig. |^ correspond to the smaller gap curve, again with 
a gratifying match between the theory and the experi- 
ment. The on-site energies have been assumed all zero 
as to implement the presence of counterions on the neg- 
atively charged backbone; the consequent induced dipole 
yet supports our fish bone construction. As a further test 
for the applicability of the present model, we have ana- 
lyzed the position of the molecular levels by comparing 
the theoretical and experimental normalized differential 
conductance, finding a fairly good accord (lower inset in 
Fig. |^). Note that the gap in the blue curve does not show 
a pronounced voltage asymmetry as observed in other ex- 
perimental curves, such as the red data. The latter mea- 
surement was performed at 3.6 K and exhibits a smaller 
gap of 0.8 eV. Here, we coped with the asymmetry by as- 
suming different voltage drops at the molecule-electrode 
junctions (bias shift of 0.06 eV) as in Ref. 

Let us now briefly discuss other possible gap opening 
mechanisms. Electron correlation may be a source for 
a gap in quasi-onedimensional systems, but for the ex- 
periment at hand band insulator mechanisms prevail on 
charge Mott one p8fl . Here, we have deliberately avoided 
the weak coupling regime since estimations of the device 
capacitance would lead to eventual Coulomb blockade 
gaps of only fractions of the observed gaps; this suggests 
that if the GC-device is in the strong coupling regime 
most of the gap is due to the molecule itself (HOMO- 
LUMO gap). Other possible gap opening mechanisms 
have been excluded from our model because of their 
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marginality to short (10 nm long) or low-tempcrature 
GC-devices. In fact (a) localization effects due to the 
sequence variability pT3| , (b) twiston motions pfl], (c) 
possible static disorder or local defects and (d) de- 
phasing p9[ are all potential causes which may concur in 
determining the absence of current in some of the recent 
experiments § @, performed at room-temperature on 
long DNA wires. 

Finally, we would also like to comment on the ob- 
served gap-width variability, even within the same sam- 
ple at different measurement sweeps. A structural fluc- 
tuation in the nucleoside distribution along the double 
helix J?!) may interfere with the it stack |}2) leading to a 
re-calibration of the overlap integrals which indeed drives 
the gap- width and induces a sharp change of the I-V pro- 
file. Moreover, the measurement process itself may in- 
duce structural rearrangement of the double helix. The 
strong electric field associated with the high nonlinear 
voltage drops can be responsible for a polarization of the 
molecule. A possible insertion of ions may result in differ- 
ent distributions of the on-site energies and the hopping 
integral, possibly locally, and leave a signature in the 
measured gap variability. Our fits show that a change 
of t± and/or tu at one site along the chain is sufficient 
to induce such a current-gap width change or a "switch" 
in the shape of the I-V curve [0], in agreement with the 
structural fluctuation hypothesis. 

Further, joint theoretical and experimental work on 
DNA molecules would definitely contribute to a better 
discrimination among the possible concomitant conduc- 
tance mechanisms, to check the influence of lead con- 
tacts on the device characteristics, and to qualify both 
the samplc-to-sample gap variability and its temperature 
dependence. New experiments may eventually be used 
to test whether low energy states can be added to the 
transmission. This could be a way to check our predic- 
tion via, e.g., changing the inter-base coupling by doping 
the molecule with metal ions || [}3[ . 

In summary, we have considered charge transport 
through a short poly(G)-poly(C) DNA molecule attached 
to nanoelectrodes by considering the hybridization of the 
7r stack with backbone states. In doing so, we have 
reached a quantitative agreement with data taken from 
the experiment reported in Rcf. |[ 
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